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Abstract 

We present the Hamiltonian analysis of the theory of gravity based on a Lagrangian density 
containing Hilbert-Palatini term along with three topological densities, Nieh-Yan, Pontryagin and 
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resulting canonical theory develops a dependence on three parameters which are coefficients of 
these terms. In the time gauge, we obtain a real SU{2) gauge theoretic description with a set of 
seven first class constraints corresponding to three SU(2) rotations, three spatial diffeomorphism 
and one to evolution in a timelike direction. Inverse of the coefficient of Nieh-Yan term, identified 
as Barbero-Immirzi parameter, acts as the coupling constant of the gauge theory. 
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I. INTRODUCTION 



Addition of total divergence terms to the Lagrangian density does not change the classi- 
cal dynamics described by it; the Euler-Lagrange equations of motion are unaltered. In the 
Hamiltonian formulation, these total divergences reflect themselves as canonical transfor- 
mations, resulting in the change of the phase space. This changes the symplectic structure 
and Hamiltonian of the system, yet the Hamilton's equations of motion remain equivalent 
to the Euler-Lagrange equations of the Lagrangian formulation. 

While the classical dynamics is not sensitive to the total divergence terms in the La- 
grangian density, the quantum theory may depend on these. The canonical transformation 
of classical Hamiltonian formulation are implemented in the quantum theory through uni- 
tary operators on the phase space and the states. However, there are special situations where 
we find topological obstructions in such a unitary implementation. In such cases these total 
divergences do affect the quantum dynamics. Therefore, to have non-trivial implications in 
the quantum theory, the total divergence terms have to be topological densities. This is a 
necessary requirement, but not sufficient. 

There are several known examples of topological terms which have serious import in the 
quantum theory. A well-known case is the Sine-Gordon quantum mechanical model [1] where 
an appropriate effective topological term can be added to the Lagrangian density to reflect 
the non-perturbative properties of the quantum theory. In this model, we have a periodic 
potential with infinitely many degenerate classical ground states. With each of these, we 
associate a perturbative vacuum state labeled by an integer n related to the winding number 
of homotopy maps S* 1 — > S 1 characterized by the homotopy group n i (5' 1 ) which is the set 
of integers Z. The physical quantum vacuum state, so called 9-vacuum, is non-perturbative 
in nature and is given by a linear superposition of these perturbative vacua with weights 
given by phases exp(m#) where angular variable 8, properly normalized, is the coefficient of 
the effective topological density term in the Lagrangian density. The physical quantities in 
the quantum theory depend on this parameter. For example, the quantum vacuum energy, 
besides the usual zero-point energy, has a contribution due to quantum tunneling processes 
between various perturbative vacua, which depends on 6. 

In field theory, we have an example of such a topological parameter 9 in the theory of 
strong interactions, namely QCD. Here also, we have infinitely many degenerate classical 
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ground states labeled by integers n associated with the winding numbers of homotopy maps 
S 3 — > S 3 characterized by the homotopy group n 3 (S77(3)) = Z. The quantum vacuum (9- 
vacuum) is a linear superposition of the perturbative vacua associated with these classical 
ground states. Associated effective topological term in the Lagrangian density is Pontryagin 
density of 577(3) gauge theory with coefficient 9. This leads to 9 dependent CP violating 
contributions to various physical quantities. However, there are stringent phenomenological 
constraints on the value of 9. For example, from possible CP- violating contribution to the 
electric-dipole moment of the neutron, this parameter is constrained by experimental results 
to be less than 1CT 10 radians. 

In gravity theory in 3+1 dimensions, there are three possible topological terms that can 
be added to the Lagrangian density. Two of these, the Nieh-Yan and Pontryagin densities, 
are P and T odd, and the third, Euler density, is P and T even. Associated with these are 
three topological parameters. In order to understand their possible import in the quantum 
theory, it is important to set up a classical Hamiltonian formulation of the theory containing 
all these terms in the action. In ref. [2], such an analysis has been presented for a theory 
based on Lagrangian density containing the standard Hilbert-Palatini term and the Nieh- 
Yan density [3]. The resulting theory, in time gauge, has been shown to correspond to 
the well-known canonical gauge theoretic formulation of gravity based on Sen-Ashtekar- 
Barbero-Immirzi real SU(2) gauge fields [4]. Here inverse of the coefficient of Nieh-Yan 
term is identified with the Barbero-Immirzi parameter 7. Thus the analysis of ref. [2] has 
provided a clear topological interpretation for 7, realizing a suggestion made earlier in [5] 
that this parameter should have a topological origin. 

The framework of [2] involving Nieh-Yan density supersedes the earlier formulation of 
Hoist [6]. Detailed Hamiltonian analysis of the theory with Hoist term for pure gravity is 
provided in ref. [7] and that including spin 1/2 fermions in ref. [8]. This discussion has also 
been extended to supergravity theories [9]. Since Hoist term is not topological, inclusion of 
matter necessitates matter dependent modification of the Hoist term so that original equa- 
tions of motion stay unaltered. On the other hand, the analysis containing Nieh-Yan density 
[2], besides explaining the topological origin of the Barbero-Immirzi parameter, provides a 
universal prescription for inclusion of arbitrary matter without any need for further modifi- 
cations of the topological Nieh-Yan term which is given in terms of the geometric quantities 
only. As elucidations of these facts, this analysis has been extended to the theory including 
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Dirac fermions in ref.[2] and to supergravity theories in ref.[10]. 

In a quantum framework, the implications of a topological term in the Lagrangian can also 
be understood through a rescaling of the wave functional by a topologically non-trivial phase 
factor. This procedure has been used for QCD [11] where, as mentioned above, the properties 
of the non-perturbative vacuum are effectively represented by a SU (3) Pontryagin density 
term in the Lagrangian. The rescaling of wave functional is provided by the exponential of 
£77(3) Chern-Simons three-form with id as its coefficient. This framework can be extended 
to the gravity theory where we have a corresponding wave functional scaling associated 
with the Nieh-Yan density. However, for the pure gravity (without any matter couplings), 
the standard Dirac quantization, where the second class constraints are implemented before 
quantization, is not appropriate. This is so because second class constraints of pure gravity 
imply vanishing of the torsion, which results in making the rescaling trivial. Instead, as 
discussed in [12], the Gupta-Bleuler and coherent state quantization methods are well suited 
for the purpose. These methods are quite general and can be used for gravity theory with 
or without matter. However, for matter-couplings leading to non- vanishing torsion, e.g. 
Dirac fermions, the Dirac quantization, as has been discussed earlier in ref. [13], can also be 
adopted for this purpose. 

Hamiltonian analysis of the first order (anti-) self-dual Lagrangian density for gravity 
including the Pontryagin density of complex SU(2) (anti-) self-dual gauge fields has first 
been reported by Montesinos in [14]. In the time gauge, the Sen-Ashtekar complex SU(2) 
connection stays unchanged, but its conjugate momentum field gets modified by the presence 
of the Pontryagin term. Recently, in [15], this analysis has also been done for gravity theories 
containing Hoist, Nieh-Yan, Euler and Pontryagin terms. This study concludes that, in the 
time gauge, real SU (2) gauge theoretical formulation is possible only if the Pontryagin 
and Euler terms are absent; the Pontryagin density can be added consistently only in the 
complex SU (2) gauge formulation leading to a canonical analysis in accordance with results 
of Montesinos [14] . 

In the following, we present a classical Hamiltonian analysis for theory of gravity based 
on Hilbert-Palatini Lagrangian supplemented with all the three possible topological terms in 
(1 +3) dimensions, namely, Nieh-Yan, Pontryagin and Euler classes. Unlike [15], in view 
of results of [2] and the remarks already made above, we shall not add the Hoist term, 
which is not a topological density. We demonstrate that, in the time gauge, we do have 
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a real SU (2) gauge theory with its coupling given by inverse of the coefficient of Nieh-Yan 
term. The canonical theory also depends on two additional arbitrary parameters, the coeffi- 
cients of Pontryagin and Euler terms in the Lagrangian density. These parameters are not 
subjected to any restrictions. A formulation of the theory presented involves the standard 
Sen-Ashtekar-Barbero-Immirzi real SU (2) connections A l a , which depend only on the coeffi- 
cient of the Nieh-Yan term, as the canonical fields. Associated conjugate momentum fields, 
instead of being densitized triads of the standard canonical theory, are modified and depend 
on the coefficients of the Nieh-Yan, Pontryagin and Euler terms. There are second class con- 
straints in the description, essentially reflecting the fact that the extrinsic curvature is not 
independent. Correspondingly, for this constrained Hamiltonian system, the Dirac brackets 
analysis is developed. Dirac brackets of the phase variables do not exhibit the same algebraic 
structure as those of the standard canonical theory of gauge fields A l a and densitized triads 
Ef; the new variables are not related to them by a canonical transformation. However, it 
is possible to construct another set of phase variables which are canonical transforms of the 
standard variables (A l a , Ef). In this framework, both new gauge fields and their conjugate 
momentum fields are modified and develop dependences on all three topological parame- 
ters. The canonical formulation described in terms these new phase variables is presented 
in detail. 

II. TOPOLOGICAL COUPLING CONSTANTS IN GRAVITY 

We set up the standard theory of pure {i.e., no matter couplings) gravity in terms of 
the 24 50(1,3) gauge connections uj/ and 16 tetrad fields as the independent fields 
described by Hilbert-Palatini (HP) Lagrangian density: 

where 

1 1 

e = det(e£) , = - efre^ = - (e^e u j - e%) , 

= { ,uj u \ j + u { l K u v]K J (2) 

and is the inverse of the tetrad field, ef, = 8% , ej = 8 T j. 
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Modifications of the gravity Lagrangian density by terms which are quadratic in curvature 
and particularly also include torsion, without altering the field equations, have a long history, 
see for example [16]. 

In (1 + 3) dimensions, there are three possible topological terms that can be added to the 
HP Lagrangian density (1). These are: 

(i) Nieh-Yan class:^] 

I NY = eEfiR^'iu) + er^D^u)e Iv D a ( U y fi (3) 

where the dual in the internal space is defined as: 

yU _ 1 jjkl v 
A = 2 KL 

and the 5*0(1, 3) covariant derivative is: D^(u)el = d e^ + 1 jC 3 . 

This topological density involves torsion. It can be explicitly written as a total divergence 

as: 

I NY = d, e[ D a {uj)e ip ] (4) 

In the Euclidean theory, as discussed in [17], this topological density, properly normalized, 
characterizes the winding numbers given by three integers associated with the homotopy 
groups n 3 (5'0(5)) = Z and n 3 (S'0(4)) = (Z, Z). 

(ii) Pontryagin class: 

I P = e^R^jHR^iu) (5) 

This is the same topological density as in the case of QCD except that the gauge group here 
is S0(l,3) instead of SU(3). Again, it is a total divergence, given in terms of the S0(l,3) 
Chern-Simons three-form: 

I P = Ad, 

For the Euclidean theory, this topological density, properly normalized, characterizes the 
winding numbers given by two integers corresponding to the homotopy group n 3 (5'0(4)) = 
(Z, Z). 

(Hi) Euler class: 

I E = e^R, uIJ {uj)R ap IJ {uj) (7) 
6 
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which again is a total divergence which can be explicitly written as: 
Ie = 49, 

For the Euclidean theory, integral of this topological density, properly normalized, over a 
compact four-manifold is an alternating sum of Betti numbers b — b 1 + b 2 — b 3 , characterizing 
the manifold. 

Now we may construct the most general Lagrangian density by adding these topological 
terms (3), (5) and (7), with the coefficients rj, 9 and respectively, to the Hilbert-Palatini 
Lagrangian density (1). Since all the topological terms are total divergences, the classical 
equations of motion are independent of the parameters 77, 8 and <fi. However, the Hamiltonian 
formulation and the symplectic structure do see these parameters. Yet, classical dynamics 
are independent of them. But, quantum theory may depend on them. 

All these topological terms in the action are functionals of local geometric quantities, yet 
they represent only the topological properties of the four-manifolds. These do not change 
under continous deformations of the four-manifold geometry. 

Notice that, while the Nieh-Yan I NY and Pontryagin I p densities are P and T violating, 
the Euler density I E is not. So in a quantum theory of gravity including these terms, 
besides the Newton's coupling constant, we can have three additional dimensionless coupling 
constants, two P and T violating (77, 9) and one P and T preserving ((f)). 

III. HAMILTONIAN FORMULATION OF GRAVITY WITH NIEH-YAN, PON- 
TRYAGIN AND EULER DENSITIES 

Here we shall carry out the Hamiltonian analysis for the most general Lagrangian density 
containing all three topological terms besides the Hilbert-Palatini term: 

£ = \ e Eft + \ I ny + {h + j Ie (9) 

where the Nieh-Yan I NY -> Pontryagin I P and Euler I E densities are given by (3), (5) and (7) 
respectively. 



(8) 
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We shall use the following parametrization for tetrad fields : 



e\ = NM 1 + N a V', el = V* ; 
MjVl = 0, MiM 1 = -1 (10) 

with N and N a as the lapse and shift fields. The inverse tetrads are: 

61 ~ N ' 6/ " Vl + N ' 
MV; = , V'V b = 5 b a , VjV? = S'j + M'Mj (11) 

The internal space metric is rj IJ = dia(—\, 1, 1, 1). The three-space metric is q a b = VjV bI 
with g = det(q ab ) which leads to e e cfet(e^) = N^Jq. The inverse three-space metric is 
q ab = yaybi^ q ab q ^ = $a Two useful identities are: 

2eE* = - V? M [7 VJi , eEfj = ^E^E^ + e A^E?] (12) 

In this parametrization, we have, instead of the 16 tetrad components e 1 , the following 
16 fields: 9 V? (M'Vf = 0), 3 M 1 (M / M / = -1) and 4 lapse and shift vector fields N, N a . 
From these, instead of the variables V? and M 1 , we define a convenient set of 12 variables, 
as: 

3? = 2eE* = e (e*e? - e'eg) = - y/q M [0 V$ , Xi = -M,/M° (13) 
which further imply: 

2eEg = - V? M [4 ^ = -Eftxn (14) 

Now, using the parametrization (10, 11) for the tetrads, and the second identity in (12), 
we expand the various terms to write: 

\ * ^Tj R,u ' » + \Iny = eTfrdvP" + - NH - N a H a - \u\ J G I3 (15) 

where we have dropped the total space derivative terms. Here tf = t]e abc D b {ui)V c i with 
e abc = e tabc and, for any internal space antisymmetric tensor, x'f] = X u + r]X IJ = X u + 



1 This parametrization differs from the one used earlier in [2] . To obtain the present parametrization replace 
eN by TV 2 in the earlier parametrization. 
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¥ijklX KL - Father, 

H = ^Z? K X%r, KL #jM = ^f K ^j b L V KL RT J ^) ~ M'DMn 
H a = eE?ji& J M = eSf^if J (c) - V*D b (u)t b 

G u = - 2D a (u;) {eXfj} = -2D a (u){e^ ta )-t a [I V J]a (16) 
where we have used the following identities: 

M'DMV = ^i a K^j b LV KL R ab IJ ^) 

V>D b {u)t] = rjeEfjRj^oo) , t a {I V J]a = -2 V D a (u) [etfj] 

Next notice that, dropping the total space derivative terms and using the Bianchi identity, 
e abc D a (u)R bcIJ = 0, we can write 

\l P + \ Ie = eh dpto" (17) 

where e a jj are given by 

(l + ifle?j = e abc {{6 + r 1 <t>)R bcI j{u) + (c/> - rj9) R bcIJ (u)} (18) 

Thus, collecting terms from (15) and (17), full Lagrangian density (9) assumes the fol- 
lowing form: 

C = nhfaP" + t a AVa - NH - N a H a - ^Gjj (19) 

with 

7r?j = e £* + e a u (20) 

In this Lagrangian density, the fields u^ IJ and ttjj form canonical pairs. Then, H, H a and 
Gjj of (16) can be expressed in terms of these fields as: 

G u = -2D a {uj)-n a I f-tl 1 V J]a (21) 
Ha = rtjRTiu) - V*D b {u)t b (22) 
H = ±- - e^) (tt$> - e g>) ^ a /V) " ^.Ht? (23) 

where we have used the relations: -D a (w)e"j = and -D a (co>)e"j = which result from the 
Bianchi identity e abc D a (u)R bcIJ (uj) = 0, and also used e b jjR ab IJ (u) = and e b u R ab IJ (u) = 
which follow from the fact that 2q (6 2 + <p 2 ) R abIJ = e abc {(9 + rj(j)) ejj -((f)- r)6) ejj}. 



Now, in order to unravel the SU(2) gauge theoretic framework for the Hamiltonian for- 
mulation, from the 24 5*0(1, 3) gauge fields ', we define, in addition to 6 field variables 
u( J , the following suitable set of 18 field variables: 

A\ ee W M« = + V u° a \ K\ ee w « (24) 

The fields A\ transform as the connection and the extrinsic curvature K l a as adjoint repre- 
sentations under the SU (2) gauge transformations. In terms of these, it is straight forward 
to check that: 

tfjdtvjp" = + <Ma )i3 = Etd t A\ + F?d t K (25) 

with 

K = - | Tfol^ = ~ \ Wi-V® = Et-^e^\A,K) + ^e^E^ Xk (26) 
Ft = 2^ + £)< = (^ + ^j{-e^E^Xk + 2e a 0i (A,K)} (27) 

where and = |e l3 ' fc e° fc as defined in (18) and effl = — rje^ are written as functions 
of the gauge field A l a and the extrinsic curvature K\ using 

RJ>) = D [a {A)K l b] - 2 - *» k K>K> 

RJM = ~\ * ijkF a b (A) + \ <P k D [a {A)K% ~ (^) KtaKi (28) 
with the SU (2) field strength and covariant derivative respectively as: 

Fi b {A) ee ^4 + 1 6^ fe ^4 , D a (A)K l b ee d a K\ + l - ^ k AiK k b (29) 

Now, using (25), the Lagrangian density (19) can be written as: 

C = E?d t Al + F t a d t K l a + t a AV^ - NH - N a H a - X - lo^Gjj (30) 

Thus, we have the canonically conjugate pairs (A l a , Ef), (K l a , F") and (V^, tf). We may 
write Gjj, H a and H of (21)-(23) in terms of these fields. For example, from (21): 

GT = \ e ijk G jk = rjD a (A)Ef + e ijk (^F fc a - t°V k ) (31) 
G b r st = G 0i = —D a (A) (E* + F^ + e^K{ { L + ±) E a k + l -F^ - t a [0 V l]a (32) 
= -D a (A)Ff + ^ k K{ | ^ + ±) El + 2 -F^ - l -e^V ak - t\ Q V i]a - i G[ ot 



10 



where the covariant derivatives are: D a {A)E\ = d a E b + rf 1 e ijk A{E b k and D a (A)F b = d a F b 
_l_ ^-i e ljk AlF^. Next, for the generators of spatial diffeomorphisms H a from (22): 

H a = E b F^ b (A) + F^D^AjK^ - K i a D b {A)F b + t b Dy a {A)V b \ - V^D b (A)t\ 

+ t b Q d [a V b ] - v a °d b t b - I (GT* + 17GJ 00 *) 

= ^d [a Ai } - A\d b E b + j^ffj - K%F b + tfa^ - Vld b t\ 

+ t b A a V b ] - V a °d b t b o + -GTK - l ~ {Gl ot + vG? ost ) K (33) 

where we have used —V*D b (u>)t b = — V^c^ + ijC^V^ + t b I V Jb uj I a J . Similarly we can express 
77 of (23) in terms of these fields. 

Now, notice that all the fields (A l a , Ef), (7Q, F") and (V^ , tf) in the Lagrangian density 
(30) are not independent. Of these, the fields V/ and t a j are given in terms of others as: 
Vl = v ! a and tf = rf with 

' K, v a = -4=Ei Xi (34) 



where is inverse of Ef, i.e., 7^£f = 5 b a , E^E* = 5) and E = det{E l a ) = ^(M )* 2 and 
T « = V e abc D b (uy c = e abc { V D b (Ay c - e^ k K{v k c + K*v°) , 

In addition, the fields F®, which are conjugate to the extrinsic curvature K l a , are also not 
independent; these are given in terms of other fields by (27) . 

In the Lagrangian density (30), there are no velocity terms associated with £0(1,3) 
gauge fields u[ J , shift vector field N a and lapse field N. Hence these fields are Lagrange 
multipliers. Associated with these are as many constraints: G u ~ 0, H a m 0, and 
H where the weak equality ~ is in the sense of Dirac theory of constrained Hamiltonian 
systems. Here from the form of G r ° l = \ e ijk G jk in (31), it is clear that these generate SU(2) 
rotations on various fields. The boost transformations are generated by G boost = G 0i , spatial 
diffeomorphisms by H a and H ?a is the Hamiltonian constraint. This, thus can already 
be viewed, without fixing the boost degrees of freedom and without solving the second class 
constraints (34) and (35), as a SU(2) gauge theoretic framework. Here, besides the three 
577(2) generators G^ ot , we have seven constraints, G boost , H a and 77. We may, however, fix 
the boost gauge invariance by choosing a time gauge. Then we are left with only the SU (2) 
gauge invariance besides the diffeomorphism H a and Hamiltonian H constraints. This we 
do in the next section. 
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IV. TIME GAUGE 



We work in the time (boost) gauge by choosing the gauge condition x% — which then 
implies for the tetrad components e° a = V® = 0. Correspondingly the boost generators (32) 
are also set equal to zero strongly G\ oost = 0. In this gauge, the Lagrangian density (30) 
takes the simple form: 

C = E«d t A\ + FfdtKi + t«d t V: - H (36) 

with the Hamiltonian density as: 

H = NH + N a H a + i e^ k ui j G r k ot + £ (V} - <) 

+ 0a (*? " Tf ) + A* - 2 ^ + ±) gg. (A, X)} (37) 

where all the fields involved are not independent. In particular, the fields V*, £° and 
depend on other fields. This fact is reflected in % above through terms with Lagrange 
multiplier fields <\> l a and \ % a . Now, in this time gauge, expressions for G[ oi , H a and H are: 

G? = V D a (A)E? + e tjk (KlF«-t*V a k ) 

H a ee E'FUA) + F»D [a (A)Kl } - K i a D b (A)F^ + t\D [a {A)V b \ - V^D b {A)t\ - ^G^K 

= ^4 - AJhM + FfW ~ K ^ P i + - V a d ^i + K ~ K a) 

H = ^L^ k E«E) {F k ab {A) - (1 + V 2 ) (D [a (A)K k - r]' 1 e kmn K^K^) } 

+ Ktf - V d a (/EG?E*) (38) 

where D a (A) is the SU(2) gauge covariant derivative. In the last line, we have used the 
time-gauge identity: t% — r o = V^EG r k ot E k . Also Ef are functions of Ef, A\ and K\: 

El = Ef(E,A,K) ee El + ?e; ( '»(A,if) (39) 

Associated with the Lagrange multiplier fields u\ 3 , A a and A in (37), we have the con- 
straints: 

G\ ot « 0, # a « , if « (40) 
In addition, corresponding to Lagrange multiplier fields and (f> % a , we have more constraints: 
VI - vi(E) « , - rf(A,A,£) « (41) 
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where, from (34) and (35), in the time gauge: 
1 



„ El, if = V e abc D b (uy c = e abc { V D b (AK - (42) 

Similarly, from the last term in (37), there are the additional constraints: 

xt = F? ~ + « (43) 

Here egj and eg, of (18), with the help of Eqn.(28), are written as functions of the gauge 
fields A % a , extrinsic curvature K l a and the topological parameters 9, besides 77 as follows: 

rf (1 + V 2 ) e a 0i (A, K) = -e abc { V (0 - rfi) - 2 V ((l - rf) - 2rfi) D b (A)K l c 

-(r, (3 - t? 2 ) £ + (3t? 2 - l) 0) e^ k KlK k c ) 

V 2 (1 + r? 2 ) eg. (A, K) = -e abc { V (9 + r?0) F^A) - 2 V ((l - V 2 ) 9 + 2#) £> 6 (A)tfj 

- ((3^ 2 - 1) - 77 (3 - V 2 ) 0) e^tf*} (44) 

From these we can construct for = + rje^ and = egj — rje^: 

a(r t ) 



e 0i 



~ ^ {^M) - (0 - rfi) D [b (A)Kl } - ^ + 2 ^ ^tf*' 

eg" = -1 e^ c {eF&A) -(9 + #) £> [6 (A)/^ - (^^) e^X*} (45) 

The constraints (43) are of particular interest. To study their effect, we note that 
(A l a , E b ) and (K l a , F b ) are canonically conjugate pairs. They have accordingly the standard 
Poisson brackets. From these, using the relation (39) expressing Ef in terms of Ef, A % a and 
K l a) as indicated in the Appendix, the following Poisson brackets can be calculated with 
respect to phase variables (A l a ,Ef) and {K l a ,F^)\ 

W*)> E b {y)] = [Ai(x), & j (yj\=5 i j %6W(x,y), 
[Ki(x), E b (y)} = , [E?(x), E b (y)} = 

These then imply the Poisson bracket relations: 

[*?(*), Kiv)] =0, Ki{y)] = -5)8 a h ^\x,y), 

[*?(*), E b (y)} = , [ X ?(x), x){y)] = (46) 

Using these, we notice that the Poisson brackets of Hamiltonian constraint H and xt are 
non-zero. Requiring [x1(x),H(y)\ ps leads us to the secondary constraints as: 

tf - (^jf) {ve ijk D b (A) (^E«E b ) + VEE^Kl} « 
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which can be rewritten as: 

f i - (^T-) ^ i^ A >c ~ e ijk Kiv k c } * 

Next, since from (41) and (42), tf w rf = e a6c { r)D b (A)vi - e ijk K J b v k }, this implies £° « 0. 
Thus we have the constraints: 

e abc {r]D b (A)vi-e ilk K J b v k c } « 

These can be solved for the extrinsic curvature K l a and recast as the following secondary 
constraints: 

i>\ = K-<(AE) « o, 

<(A,£) ee ^ e ijh E J a D b (A)El 

~ ^y cd {E k b D c (A)E d + ElD c {A)E k d - V k E™D c {A)E™} (47) 

These are additional constraints and have the important property that these form second 
class pairs with the constraints xt °f (43): 

[*?(*), $'(y)] = -SSfil6®(x,y) (48) 

To implement these second class constraints, and V'a) we nee d to go over from Poisson 
brackets to the corresponding Dirac brackets and then impose the constraints strongly, ip a = 
(which also implies t% = 0) and xt — 0, in accordance with Dirac theory of constrained 
Hamiltonian systems. As outlined in the Appendix, the Dirac brackets of fields A\ and Ef 
turn out to be the same as their Poisson brackets; these are displayed in (A. 20). On the 
other hand, those for (A a , E®; K a , F?) are different; these have been listed in (A. 22) and 
(A. 23). 

Finally, after implementing these second class constraints, we have the Lagrangian density 
in the time-gauge as: 

C = E?d t A a + F t a d t K a - H (49) 

with the Hamiltonian density 

H = NH + N a H a + - e ijk u{ k G r k ot (50) 
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and a set of seven first class constraints: 



GT = V D a {A)E? + t l3k K{F a k « 

H a = E^JA) + F'D^A)^ - K i a D b (A)F b - rj' 1 G^K^ « 
H s ^ C «*^F* (A) - f i±f 



EEfE^m + -d a (VEGfEZ) « (51) 



with £7? in the last equation given by: Ef = Ef(E, A, K) = Ef + 2 e^ v) (A, K). The fields 
(A l a , Ef, K l a , F- a ) have non-trivial Dirac brackets as listed in (A. 22) and (A. 23). The second 
class constraints xt an d V'o are now se ^ strongly equal to zero: 

Kl = kI(A,E) = \ e^ k EiD b {A)El - ^ E k a e bcd {E k b D c {A)E^ + ElD c {A)E k d 



V k E™D c {A)E™} 



Ft = 2[ V + -)e a 0i (A,K) 



(52) 



In writing the Hamiltonian constraint H in (51) from (38), we have used the identity: 



(53) 



which holds due to the time gauge relation EEfG™ 1 
K l a = K l a (A,E) imposed strongly. 

To evaluate the effect of generators (51) on various fields, we need to use the Dirac 
brackets instead of the Poisson brackets. For example, for the SU(2) gauge generators, 
using the results listed in the Appendix, we obtain: 



e abc E l b K l c with the constraints 



D 



^ k E a k S^(x,y) , 



[GT*(x), Ai(y)] D = - V (5^d a +r ] - 1 e ikj A k ) S^(x,y) 



(54) 



reflecting the fact G r ° l are generators of SU(2) transformations: A\ transform as the SU{2) 
connection and fields Ef as adjoint representations. Besides, the fields Ff, K\ and Ef also 
behave as adjoint representations under SU{2) rotations: 



D 



G7\x\ F«{y) 
[G?(x), Ki(y)] D 



e Hk K k sW( X ,y) 

e ijk E a sW( x ,y) 



(55) 
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Similar discussion is valid for the spatial diffeomorphism generators H a . The Dirac brack- 
ets of H a with various fields yield the Lie derivatives of these fields respectively, modulo 
SU(2) gauge transformations. 

As stated earlier and demonstrated in the Appendix, Dirac brackets for the fields 
(A* , Ef\ K l a , F t a ) are different from their Poisson brackets (see (A.21), (A.22) and (A.23)). 
This is so because the transition from Poisson brackets to Dirac brackets, except for some 
special cases, in general, does not preserve canonical structure of the algebra [18]. When the 
second class constraints are imposed strongly, the algebraic structure of the Dirac brackets 
of phase variables (A l a , Ef) of the final theory is different from those of the phase variables 
(A l a , Ef) of the standard canonical theory. Thus the variables (A l a , Ef) are not related to 
(A l a , Ef) through a canonical transformation. However, it is possible to construct a set of 
new phase space field variables whose Dirac bracket algebra has the same structure as that 
of the standard canonical variables (A l a) Ef). 

In fact, in general, for theories with second class constraints as is the case here, instead of 
the ordinary canonical transformations, what is relevant are the Gitman D-transformations, 
which preserve the form invariance of Dirac brackets and equations of motion [19]. Thus, 
in the present context also, new phase variables can be constructed through these D- 
transformations. These transformations change both the gauge fields as well as their conju- 
gate momentum fields. This procedure finally leads to the phase variables: 

oo 1 oo 1 

lb* ft' 

n=0 n=0 

where 

D'^ a (x) ee Ef(x) , CT ] \x) ee A\{x) (57) 

and other D'M and C'^ are recursively constructed using Dirac brackets as follows: 

D't +1)a {x) = J d 3 z F*(z) 

C'^\x) = Jd 3 zF?(z)[K l b (z), Cf*(s)] D - 2 - J dhe% )a (z)[A l b (z), C**{x)] D 

n = 0, 1, 2, 3, (58) 

In particular, 

D?\x) = Ft(x) - 2 - 4f (A, K; x), Cf(.) = F*(z) K ^ Z) (59) 
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K l b (z), D't )a (x) 



d 3 ze^ a (z) 



J n 



The new variables (A 1 ^ Ef) are functions of the phase variables {A\, Ef, K l a , Ff) of 
the theory described above and can be checked to satisfy the Dirac bracket relations: 

[K(x), Ef(y)] D = 5{5 b a 5^(x,y), [A«{x\ A^y)] D = 0, [Ef (x), E?(y)] D = (60) 

As is expected under D-transformations, these relations reflect the fact that the algebraic 
structure of Dirac brackets for the fields (A % a , Ef) as represented by (A. 20) has been pre- 
served. After the second class constraints, xt an d ip l a , are implemented, (A'*, Ef) are related 
to the phase variables (A l a , Ef) through an ordinary canonical transformation. 

We have not presented many details of the construction of these new phase variables 
above. Instead, in the next section, we shall present, through an equivalent procedure, an 
elaborate construction of the new phase variables in the theory where second class constraints 
are already imposed strongly. This will be done by a direct canonical transformation of the 
phase variables (A l a) Ef) of the standard canonical theory. The new canonical variables so 
obtained will be shown to be equal to the fields A'* and Ef above, when the second class 
constraints xt an d Va are imposed. 

V. CANONICAL TRANSFORMATIONS AND NEW PHASE VARIABLES 

Adding the Nieh-Yan term to Hilbert-Palatini Lagrangian density, in the time gauge, 
leads to a change of phase variables [2], from the ADM variables Ef) to new variables 
(A\, El). This change is just a canonical transformation. Further inclusion of the Pontray- 
gin and Euler densities results in a theory which can also be described in terms of canonically 
transformed phase variables. In the following, we shall develop such a description explicitly. 

We start with the standard canonical theory constructed from the Lagrangian density 
containing the Hilbert-Palatini term and the Nieh-Yan density as in (9) with 6 = and 
= 0. This is described, after partial gauge fixing (the time gauge), where the second class 
constraints are imposed, in terms the SU(2) gauge fields A\ and their conjugates, densitized 
triads Ef by the Lagrangian density 

n = ^ e ijk u\ j G r k ot + N a H a + NH 
Gl°\A,E) = i]D a (A)E? , H a (A,E) 



= EtFUA) - rf\\G\ 
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H(A,E) = ^ e^EfE* (F k b (A) e kmn K™K^ + ± d a (VEcrfE*) (61) 

where the extrinsic curvature K l a (A, E) is given in terms of A l a and Ef through (A. 7). Canon- 
ical pairs of the phase variables (A l a , Ef) obey the standard Poisson bracket relations: 

[Ai(x), E){y)] = b) 5 b a 6<*>(x,y) , [A\(x), A{{y)] = , [Ef(x), E){y)] = (62) 

Next, we add the Pontraygin and Euler densities (6, 8), which are total divergences, 
jl p + \l E = J M , to C 1 above. The resulting Lagrangian density, ignoring the spatial 
derivative part, is 

C 2 = E^Ai + d^-H (63) 

Inclusion of the time derivative term here is equivalent to a canonical transformation on the 
phase space which can be constructed using J*. For this purpose, we first express J* as a 
function of the phase variables A l a and Ef: 

J\A,k(A,E)) = ^e abc ^A\Fl c -^e^A\A{A k ^ 



+ ^{(l- v 2 )e + 2 V <f ) }e abc ^D [b (A)Kl ] 

+ A {(3^ - 1)9 - V (3 - V 2 )<P} ^ c ^ k ^A< (64) 

Generating functional for the canonical transformation is: 

J{A,E) = J d 3 z J t (A(z),K(A(z),E(z))) (65) 

which has functional dependence on both gauge fields A l a and their conjugates Ef. Follow- 
ing the standard procedure, J generates the canonical transformations, (A % a (x), Ef(x)) — > 
(A l a (x), £f(x))i where the new phase variables are given in terms of Poisson bracket series 
as follows: 

Al(x) = A\{x) + [j,Ai(x)} + 1 [J, [j,A\(x)]\ + i [J, [J, [J,A^x)}}} + 



2! L ' L ' ay /n 3! 
= e J Ai(x) e~ J 

£f(x) = Ef(x) + [J, Ef(x)} + 1 [J, [J, Ef(x)}} + 1 [J, [J, [J, £?(*)]]] + 

= e- 7 (66) 
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Alternately, these relations may be represented as: 

A ^ = E ^ . w = E h A (n)a ' 

with 



n=0 



ra=0 



n! 



(67) 



cl 0)l = K(x), 



D^ )a (x) = Ef(x) 



J, D [ r 1)a (x) 



, n = 1, 2, 3, 



^ n)l (*) = [.7, c^- 1 ^*)] , A (n) » = 

The various terms can be evaluated recursively through the following formulae: 
C^\x) = J dh (D^\z) [A[{z\ C^(x)]-Ci 1)l (z) [E*(z), Cf^s)]) 



(68) 



D) (x) 



d 3 z (Dl 1)b ( 



A{(z), Dt 1)a (x)] -Cl 1)l (z) [E*(z), Dt 1)a (x) 
n = 1,2,3 



(69) 



where, using J from (64, 65), 

C^\x) = [J, A\{x)] 



= [J, £?(x)] = 2C a (« ;a ;) 



77 



6E?(x) ' 



(70) 



Here e^ a (n] x) = e^re; x) + 77^ (k; x) and the argument k is to indicate that these functions 
are given by (44) with K l a replaced by K l a (A, E) = A\ + f l a {E) where f l a {E) are as in (A. 7). 

By repeated use of Jacobi identity, it can be checked that the functions and of 
(68) satisfy the Poisson bracket relations: 



E 

1=0 

n 

E 



E 



1 


l\(n- 




1 




l\(n- 


l)\ 


1 




l\(n- 


1)1 



c^(x), ct l)3 {y)} = 

D? a (x), D[ n - l %) 



= 



0, n = l,2,3, 



(71) 



The Poisson bracket relations (62) imply, by construction, same Poisson brackets for the 
new variables (66): 



£j(v)} = % % t (3) (^y) , Mv)} = o , [^(s), £?(</)] = o (72) 
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where the Poisson brackets are evaluated with respect to the phase variables (A l a , Ef). This 
can be readily checked by using the identities (71). 

For a general analytic function P(A, E) of the phase variables A l a and Ef, the following 
relation holds: 

P(A, £) = e J P(A, E)e~ J = P(A, E) + [J, P(A, E)\ + I [J, [J, P(A, E)}} 

+ ^[J,[J,[J,P(A,E)]]] + (73) 

Further J of (64 , 65) written as a functional of (A l a , Ef) and (A l a , £f) is form invariant: 

J(A,£) = J(A,E) (74) 



The converse relations expressing A\ and Ef L in terms of the transformed variables are: 



Ai(x) = A\{x) - [j,Ai(x)} + 1 [J, [j,Ai(xj\] - ± [J, [J, [J,Al(x)]]] + 



= e~ J A\{x) e J 

E?{x) = £f{x) - [J,£?(x)] + | [J , [^(s)]] - | [J , [J , [^(s)]]] + 

= e^ 7 (75) 

where J is written as a functional of ^ and £f (refer (74)) and Poisson brackets are 
evaluated with respect to these new variables. 
Next, we evaluate the following: 

/oo n „ 

d 3 x = J2 F(n) > F(n) = E 7i7773nT / ^ ^^^(a;) 

It is straight forward to check: 

F<°> = | rf 3 x ^(x)^4i(x), F« = 9 4 GW + d t J 

F (n) = - d t G {n ^\ n = 2, 3, 4, (76) 

where 

n . „ 

G W ss [J , G (n_1) ] = £ y rf 3 x A (/)a (^)^ +1 ^(^) , n = 1, 2, 3 

G {0) = | rf 3 ^ ^(x)C7W*(x) (77) 
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To obtain this result, the following helpful identities may be used: 



n-1 



E Z!( w -i-Z)! / d " X { D t l)a SC^ - C^5D? a ) = 0, n = 2, 3, 4 (78) 

which can be derived recursively by taking Poisson brackets with J . 

Further, using expression for C^ 1 from (70) and eqns. (A. 9, A. 10) of the Appendix, the 
following relation can be obtained: 

Et(x)Ci^) = (l + V 2 ) e ijk d a {e b 0l (K;x)El(x)E a k (x)) 

which in turn implies G (0) = / d 3 x Ef(x)Ci 1]i (x) = and hence all the of (77) are zero, 
thus leading to the result: 

J d 3 x £?{x)d t J¥ a {x) = f d 3 x E«(x)d t Ai(x) + d t J (79) 

Since the generating functional J~, as given by (64, 65), is invariant under small SU(2) 
gauge transformations and spatial diffeomorphisms generated respectively by G^ ot (A, E) and 
H a (A,E) of (61): 

[j, G™\A, E)] = , [J, H a (A,E)]=0 

Consequently, G\ ot and H a written in terms of the phase variable (A\, Ef) and (A\, £?) are 
form invariant: 

G\ ot {A,£) = e J G r t ot (A,E) e- J = G r t ot (A,E) 
H a (A,£) = e J H a (A, E) e' J = H a (A, E) (80) 

On the other hand, for the Hamiltonian constraint we have: 

H(A,E) = e ~ J{A£) H(A,£) e J{A ' e) 

= H(A,£) - [J,H(A,£)\ + ^[J,[J,H(A,£)\] 

- \ [ [JAJAJ,H{AE)\\\ + (81) 

where the Poisson brackets are with respect to phase variables {A l a , £f). 

This detail discussion, finally allows us to write the theory based on the Lagrangian 
density (63) in terms of the new phase variables as: 

c 2 = £«d t Ai - n 

% = ^ k 4 j G r k ot (A, 8) + N a H a (A, £) + NH(A, £) (82) 
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where 



or {AS) = r)D a (A)£? 
H a (AS) = s>r ab {A)-^G7\AS)^ a (A£) 

H(A,S) = e- J ^ H(A,£) e J ^ (83) 

The new variables (A\, £?) obtained here are related to the variables E[ a ) of eqns. 
(56-60) of Sec. IV derived by the Gitman D-transformations. When the second class con- 
straints xt an d Va there are implemented, (A'*, E'®) collapse to (A l a , £?): 

A>:(x = 0, V = 0) = Al , E>?(x = 0,^ = 0)= St (84) 

This is so because each of the terms in (56) and (67) coincide: 

C'^\x = 0^ = 0) = Cf* , D't )a {x = 0,^ = 0) = D { t )a (85) 

This completes our discussion of the canonical transformation to new variables (A l a , Sf) 
obtained by adding the Pontraygin and Euler densities to the standard canonical theory of 
gravity described in terms of the phase variables (A l a , E?). 



VI. SUMMARY AND CONCLUDING REMARKS 



We have developed the canonical Hamiltonian formulation of gravity theory with all 
the three topological terms of the Lagrangian density (9) as an SU(2) gauge theory with 
Barbero-Immirzi parameter 7 = i]^ 1 as its coupling constant. In time-gauge, the theory 
containing only the Nieh-Yan topological term {9 = 0, = 0) developed earlier in ref. [2], 
is described by real SU(2) gauge fields, A l a = oo° a l + r/uj^ 1 , and densitized triads E? as their 
conjugate momentum fields. This coincides with the standard SU(2) gauge theoretical 
canonical formulation of the theory of gravity [4]. When the Pontryagin and Euler terms 
are also included, there is a formulation of the theory which retains the gauge fields A\ 
(independent of the topological parameters 9 and 0) as the canonical fields, but, their 
conjugate momentum fields are modified from E? to Ef = Ef — 2^ _1 e^A, K) + le^A, K) 
developing dependence on 9 and 0. Further, for the case with 9 = and = 0, the 
momentum conjugate to extrinsic curvature K l a is zero. Here, in the most general case, it 
is non-zero, represented by Ff which depends on other fields through the Xi constraints 
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(43). In addition, it also depends on the topological parameters 9 and 0. Associated with 
Xi , we have a set of secondary constraints ip % a of (47) which expresses the fact that extrinsic 
curvature K l a is not an independent field. These constraints, ip l a ), form second class 

pairs which are implemented by going over to the Dirac brackets from Poisson brackets. 
The theory is described by seven first class constraints: the £77(2) gauge constraints G\ ot : , 
spatial diffeomorphism constraints H a and Hamiltonian constraint H as listed in (51). In 
this formulation, however, the Dirac brackets for the phase variables (A l a (x), Ef(x)) do not 
possess the same algebraic structure as those for the canonical variable (A l a (x), Ef(x)) of 
the standard theory. Even after the second class constraints, xt — 0> — 0> are imposed, 
there is no canonical transformation that relates the set (A l a , E?) to (A l a , E?). However, it 
is possible to construct another Hamiltonian formulation in terms of new canonical variables 
(A l a , £f) which indeed are related to the standard variables (A % a , Ef) through a canonical 
transformation. Here both the gauge fields as well as their conjugate momentum fields, as 
represented in (66), are changed and these depend on all the topological parameters, 77, 9 
and 0. 

From this classical Hamiltonian formulation described in terms of (A l a , Sf), we can go 
over to the quantum theory by replacing the Poisson brackets by commutators of correspond- 
ing operators in the usual fashion. We already have some evidence that Barbero-Immirzi 
parameter rj~ l is relevant in the quantum theory. For example, it appears in the spectrum 
of area and volume operators [20] and also in the black hole entropy [21]. How other param- 
eters, 9 and 0, will be reflected in the quantum theory is an open question requiring deeper 
study. 

The analysis presented in the present article is for pure gravity without matter couplings. 
Inclusion of matter, such as fermions, spin 1/2 or spin 3/2 (supergravity), may be achieved 
through standard minimal couplings. All the topological densities in the Lagrangian are 
described in terms of geometric quantities only. Their presence does not change the classical 
equations of motion even with matter. A Hamiltonian formulation, in the time gauge, can 
again be set up in terms of a real £77(2) gauge theory with rj~ l as its coupling constant. 
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Appendix: Poisson and Dirac brackets 



In the time-gauge Lagrangian density (49), the fields (A l a , Ef) and (K^, F?) are canonical 
pairs which have the standard Poisson bracket relations: 



5}5 b a S^(x,y) , [K(t,x), F b (t, y)] = 5)5 b a 5^(x,y) (A.l) 



and all other brackets amongst these fields are zero. Thus the Poisson bracket for any two 
arbitrary fields P and Q is given by: 

\p(x) Q (y)\ = [ dh (EKm± _ i^L mi) 

1 { h Qm J \5AKz) 5E? {z) 5E?(z) SAl(z) ) 



+ 



, 3 , SP(x) SQ(y) 5P(x) SQ(y) 
a z 1 



SK*(z) 8F?{z) 5F-{z) SKfa) 
,-i ~ a (v) i 



(A.2) 



/From these, using Ef = Ef(E,A,K) = Ef + 2r]~ l e^>{A,K), we have the Poisson 
bracket relations 

K0*0> E j(y)] = E%y)] = 5}5 b 5^(x,y), [K(x), E b (y)] =0 (A.3) 

Using the expressions for e^A, K) and e^A, K) as functions of A\ and K\ as in (44), the 
following relations obtain: 

'9 + T](f)' 



~E*(x),E b (y) 


2 
V 


[E?(x),e b f(y)_ 


4 
rf 


'_mx),E b (y) 


2 
V 







nbc 



9D l l - 



e ikj K k }S {3 \x,y) 



eof(x),E b (y) = e^(x),E b (y) 



2 abc W - ( °-^\ e^K k c \ S^(x,y) 



V 



V 



(1 + rf) [e a 0i (x), E b (y)} = (1 + rf) e a 0i (x), E b (y) 



- e 

V 



abc 



(9 + 7)^- 



rf)9 + 2rj(f)\ 
V J 



e ikj K k \s®(x,y) (A.4) 
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where the SU(2) gauge covariant derivative is: D l J = 5 tj d c + r] 1 e lkj A k . These Poisson 
bracket relations imply for Ef = Ef{E, A, K) = Ef + 2^ e^{A, K): 

[E?(x), E]{y)] = (A.5) 

Now, using these Poisson bracket relations along with (A. 3), yields: 

E){y)] = [A\{xl E){y)] = 5 h a 5) 5 3 (x,y) (A.6) 

where K a (E,A) is given by (47) and can be rewritten explicitly as: 

K a (A,E) = A a + f a (E) 



fi(E) = { e^ k Eid b E b k - ± E k a e bcd (Eld c E k d + E k d c E d - 5 ik E l b d c E l d ) 



\ e l3k Eld b E b k - ± 
= ~ V E{ e bcd (vid c v d - \ F v r b d c v^ (A.7) 
with v\ = E l a /\fE. It is straight forward to check that /* satisfy the identity: 

e abc {d b El - d b (\n VE) El - rf\^ k f b E k } = 
Equivalently, this relation can also be written as: 

d a Ef - v-WfiEl = D a (A)Ef ~ r^ l3k <Et = (A.8) 
These relations can be used to calculate the variation 5f l a to be: 

& fa = ^ b 5E b - rj d c (A ab cU 5E b ) + \ (d c A ab al ) 5E b (A.9) 

with 

Si = - (E l a r b + Elf a ) + 1 (E*fl + E[f a ) - \ E™E™ (Eff c + Eftf) 

+ \ (E™E l h E1 + E™E a Ef) f™ + (E l b E l a - E\E[ + 5 li E n a E%) E C JT 

- I e mk [d c E™E l b - d c E l b E?) E c k -\ e lmk {d c E™E a - d c E a E™) E% 

+ \ e ak {d c E™E™ - d c E™E™) E c k 

f 1 1 \ 

A cil I Alk rpm rpm Amk rpm rpl . Irak rpm rpi \ rpc ( \ ia\ 

A ab = I e E a E b - - e E a E b + - e E b E a \ E k (A.1U) 

Notice that S ab and A ab cd are respectively symmetric and antisymmetric under the inter- 
change of the pair of indices (a,i) and (b,l): 

S a \ = S l ba , A ab cU = - A ba cli (A.ll) 
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These properties, immediately, lead to the relation: 

Sfiix) _ 5fl(y) 



(A.12) 



5E\{y) 5E?(x) 

Next, using Xii x ) = E^(x) ~ ' ^oi( x ) fro m (43), equations (A.l) also imply the 
following: 

2(1 + 77 2 ) 



X a i{x),my) 



(i 



~^ + 2 ^ ) e^K^S^(x,y) 



2(1 + V 2 
V 



eUx),PKy) 



= { ((1 - ^ + 2,0) Z>« + ( W~ 1)^(3-^ e^} ^(x, y) 

[#(*), ^j(y)] =0, ^(y)] = -«^(x,y) 

(W) e b 0j (y)} = (1 + rf) [Ft(z), e b 0j (y)] 

= 2 - { ((1 - , 2 )0 - 2rfi) D» + ~ ^ + ^ ~ ^ e^K*} 5*\x, y) 
(l + V 2 ) [*?(*), e b 0j (y)} = (1+, 2 ) lifts), e a 0j (y) 



[Xi(x), e 



v e 



abc 



((1 - V y + 2,0) D« + ( W~ 1)^(3-^ e^} ^(s.y) 



= (TO, ef )b (y)] 

^ e abc {(0 + - (^ ((1 ~^ + 2# N ) e^ c fc |> ^(x,y) (A. 13) 



which further imply: 

^(y)] = , [x?(x), «j(y)] = 0, [xUx), X){y)] = (A.14) 

For ip l a = K l a — K l a (A,E) as given by (47), using (A. 3) and (A. 7), we have the following 
useful relations: 



K(x), E b (y)} = 

K(x), Ai(y)] = 

[ti(x), El(y)} = 

E{y)] = 



[<(x),E b (y)}-- 
[<(x), Ai(y)} 
[<(x), El(y)] 
[<(x), E(y)} : 



-6^ 5 3 (x,y), 

$<(x) 

SE b (y) 
= EiEl5 3 (x,y), 



EEl8 3 (x,y) 



(A.15) 
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The Poisson bracket relations among xt and V'oj obtained by using the properties listed 
above, can be summarized as: 

[rf(s),X$(y)] = , [x?(x),i>i(y)} = -S^j 6®(x,y) , [#(s),$(v)] = (A.16) 
where the last equation follows from the relation: 
R(x), <{V)] = Wx), fl(y)] + [P a (x), Ai(y)} = - ^ = (A.17) 

Here the Poisson brackets involving f l a {E) are calculated by using their expressions as func- 
tions of E l a as given by (A. 7). The identity (A.17) further implies the following Poisson 
bracket relations: 

RW. <{V)\ + [D [a {Ay b] {xl Ai(y)]=0, 
Fjx), D [c {A)K^{y)] + [D [a (A)4](x), i&(y)] = 0, 

d^aK^^^akM + [FM,z jmn K?(yK(y)} + b ikl <(M(^FL(y)} =o 



= o, 



(A.18) 



To implement the second-class constraints xt ~ and ~ 0, we need to go over to the 
corresponding Dirac brackets and then put xt — an d = strongly. From the Poisson 
bracket relations of these constraints (A.16), the Dirac bracket of any two fields C and D 
can be constructed to be: 

[C, D] D = [C, D]-[C, x] D] + [C, V] [x, I>] (A.19) 

Using the Poisson bracket relations listed above, it is straight forward to check that the 
Dirac brackets amongst A\ and Ef are the same as their Poisson brackets: 

= 0, [A\{x) ) A{{y)] D = [Ai(x),Ai(y)} = 



[E!(x),E b J (y)] D = [E«(x),E b (y 
[A i a (x),E b j (y)] D ^ [Al{x),E){y)] 

Also we note that, 

[K(x),E b j (y)] D =[Ki(x),E b j (y)] D 

[Ki(x) 1 ^(y)] D =[Ki(x) 1 Al(y)] D 



5 b a 5)5\x,y) 



(A.20) 



[<(x),E b (y)] = [Ai(x),E b (y)} = 5 b a 5) 5 3 (x,y), 
[<(x),Al(y)] = [r a (x),Ai(y)} = - 
[<(x),4(y)] 

[A\(x)Jl(y)\ + [f a (x),Al(y)]=0 (A.21) 
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where in the last terms of second and third equations, the Poisson brackets are to be evalu- 
ated using (A. 7) which express f l a (E) as functions of E\. 



The Dirac brackets of (A l a , Ef) and (Ef, E b ) are not same as their Poisson brackets: 



K(x), E b (y) 



D 



rj 



2 r 



= 8$6\x,y)-- A* a (x), ef 

D V L 



D 7] z 



+ [e i ™K™(x)K n d (x), F> f (yj\) 



4(fl 2 + 



^acd^bef 



(A.22) 



D [e (A)"b(x), D { MW f] (y) 

Here the argument k in eg^/c) and eg^/c) is to indicate that these are as in (44) with K l a 
replaced by k\ which in turn are given by (A. 7) as functions of A\ and E®. Further, here in 
the second equation, we have used: 



Also, 



K(x), F b {y) 



E?(x), y) 
2(1 + 77 2 ) 



+ 







D 



D 



V 

2(1 + rj 2 ) 
V 



V 

2(1 + rj 2 ) 
V 



[E?(x), e^K-y)] 



F?(x), F b {y) 



J D 



V 

4(1 + T] 2 
r] 2 



[e 2 + 2 ) 



2^ ^acd^bef 



D [e (A)^(x), D [e (AW f] (y) 



Et(x),F b (y) + F«{x),E h {y) 



D 



_ 8(fl 2 + 2 ) 
2(1 + ?7 2 



[eSi(«;^)> e b oj( K ;y)] 



D [c (A)^(x), D [e (AW f] (y) 



D 



7/ 



KW. = [««(*)» ( A - 23 ) 
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